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I. INTRODUCTION 

The phrase strange nonchaotic attractor Q (SNA) 
refers to an attractor which is nonchaotic in the sense 
that its orbits are not exponentially sensitive to pertur- 
bation {i.e., none of the Lyapunov exponents are posi- 
tive), but the attractor is strange in the sense that its 
phase space structure has nontrivial fractal properties. 
Past studies indicate that SNA's are typical in nonlin- 
ear dynamical systems that are quasiperiodically forced. 
Here by a typical behavior we mean that the behavior 
occurs for a positive measure set of parameter values. 
Alternatively, if parameters are chosen at random from 
an ensemble with smooth probability density, then the 
probability of choosing parameters that yield a typical 
behavior is not zero. The description of a behavior as 
typical is to be contrasted with the stronger statement 
that a behavior is robust. In particular, we say a behav- 
ior of a system is robust if it persists under sufficiently 
small perturbations; i.e., there exist a positive value 5 
such that the robust behavior occurs for all systems that 
can be obtained by perturbation of the original system 
by an amount less than 5. Thus all robust behaviors are 
also typical, but not vice versa. 

With respect to SNA's, examples where they are 
typical but not robust have been extensively studied 
IHj §7 @7 Hi- All example of this type is the quasiperi- 
odically forced circle map given by the system , 

dn+i = [dn + mod 27r, (la) 
■i^n+i = [<y5n + 1^1^ + £ sin + Ccos 0„] mod 27r, (lb) 

where 17 = Lo/iTi is irrational. Other examples of typi- 
cal nonrobust SNA's involving differential equations have 
also been studied [|[ |j. Numerical evidence H, ^ and 
analysis based on a correspondence ||, |^ with Anderson 
localization in a quasiperiodic potential leads to an un- 
derstanding of the typical but nonrobust nature of SNA's 
in these examples: In particular, it is found that SNA's 
exist on a positive Lebesgue measure Cantor set in pa- 
rameter space. In the case of Eq. (|l|), for example, con- 
sider the rotation number 

W = lim ((^„-</7o)/(27rn), (2) 



where for this limit Lpn is not computed modulo 27r. For 
fixed w, e > 0, and C > 0, a plot of W versus ijJ^ yields 
an incomplete devil's staircase, a nondecreasing graph 
consisting of intervals of w^p where W{ijJip) is constant and 
with increase of Wiuj^) occurring only on a Cantor set 
of positive measure. For small e, the values of uj^p on the 
Cantor set correspond to orbits that are three frequencies 
quasiperiodic, but for larger e they correspond to SNA's. 
Because an arbitrarily small perturbation of ui^p from a 
value in the Cantor set can result in a value of uo^p outside 
the Cantor set, these SNA's are not robust. On the other 
hand, because the Cantor set of uj^p values has positive 
Lebesgue measure ('positive length'), these attractors are 
typical for (^. 

Other studies suggest that there are situations where 
SNA's are robust |l , |, @, |, |, |l^, |l| . The experi- 
ment of Ditto et aL[q] on a quasiperiodically forced mag- 
netoelastic ribbon produced evidence of a SNA, and the 
existence of this SNA appeared to be stable to parameter 
perturbations. The original paper where the existence of 
SNA's in quasiperiodically forced systems was first dis- 
cussed 1^ gives numerical evidence of robust SNA's. In 
addition, the effect of quasiperiodic perturbations on a 
system undergoing a periodic doubling cascade has been 
investigated, and evidence has been presented suggesting 
that, after a finite number of torus doublings, a robust 
SNA results @,§. 

Thus there seems to be two types of SNA's: typical, 
non-robust SNA's, and robust SNA's. In this paper we 
study a class of models exhibiting robust SNA's. The 
model class that we study is particularly interesting be- 
cause it allows the possibility of rigorous analysis. In 
particular, we are able to prove, under the mild hypoth- 
esis that a certain Lyapunov exponent is negative, that 
the attractor is strange and nonchaotic. Since other cases 
of SNA's are likely to be accessible only to study by nu- 
merical means, it is worthwhile to investigate our, more 
well-understood models, numerically. By doing this we 
gain insight into the applicability and limitations of nu- 
merical techniques for the study of SNA's. 

In this paper we consider quasiperiodically forced maps 
which can be motivated by consideration of a system 
of ordinary differential equations in the form dx/dt = 
F(x,^, 61(1), 61(2), . . . , 0(^)), where F is 27r periodic in the 
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angles ^ and which are given by ^ = uj^t + ^q, 

9^^^ = ujg{i)t + and w^, , • • • ,we(jv) are incom- 
mensurate. Sampling the state of the system at discrete 
times tn given by ^ = 2n7r, we obtain a mapping of the 
form, 

6f+, = [0^^ +uj^'^] mod27r, (3a) 
x„+i=F(x„,0«,0(2)^... ^gW)^ (3^^) 

where x„ = x(t„), cj*-*) — 2Trujg(i) /uj^, and 
there exist no set of integers (m(°\ m'^^\ • • • ,m^-^^) 
for which E^j^m'^^^w*-*-' = 2Trm^^\ aside from 
(mW,m(i),--- ,m(^)) = (0,0,--- ,0). 

For the map (^, the simplest possible attractor is an 
TV-dimensional torus, x = i {9'-^\ 9'-^\ ■ ■ ■ , 6l(^)). In this 
paper, we consider the case where an attracting {N + 1)- 
dimensional torus exists, and the dynamics on the torus 
is given by 

0^U = [Oi:^+co^'^] mod27r, (4a) 
= + q^'^eip + + . . . + q(^Hir (4b) 

+ P(^„,C),0f mod 2^, 

where P is periodic in all its variables and q^^\ ■ ■ ■ j 
q^^'^ are integers. We are particularly interested in the 
case that (^) is invertible, so that no chaos is possible, 
and when at least one g^*^ is nonzero, which as we will 
see prevents the existence of an attracting iV-torus. 

In Sec. II we examine the simplest case where = 1 
(^(*) 9). Section II. A presents numerical experiments 
and rigorous analysis of this two-dimensional map model. 
In particular, we prove (subject to a mild hypothesis on 
the negativity of a Lyapunov exponent) that, for our 
class of maps, the information dimension of the SNA is 
one (Z?i = 1), while its box-counting dimension is two 
[Dq = 2) psj. Thus we rigorously characterize the na- 
ture of the strangeness of the SNA's for our model. [In 
a previous work [ p^ it was argued (nonrigorously) that 
Di = 1 and ZJq = 2 for the two-dimensional SNA map in- 
troduced in Q.] We conjecture that Di = 1 and Dq = 2 
typically holds for SNA 's of two-dimensional quasiperiod- 
ically forced maps {i.e., maps of the form 9n+i — {9n+uj) 
mod 2tt, ipn+i — F{ipn,9n) with Q = uj/2tt irrational). 
Also, in Sec. II. A we present numerical experiments on 
dimension calculations of Di and Dq, and of the Lya- 
punov exponent for our map. Section II. B investigates 
the dynamical origin of SNA's as a limit as fl approaches 
its irrational value through an infinite sequence of finer 
and finer rational approximations (RA's) [ p^ . It turns 
out that this technique yields substantial insight into the 
structure of SNA's, as well as additional understanding 
of why Di — I and Dq — 2 applies. 

Section III considers higher dimensional maps. In par- 
ticular. Sec. III. A considers the case where x in Eq. (||b) 
is two-dimensional and = 1, while Sec. III.B considers 

> 1 with x a scalar angle variable (as in Sec. II). For 
the map of Sec. III. A, we consider one component of x 



to be an angle variable and the other component to be 
analogous to a radial variable. Thus, if, on the attractor, 
the radial coordinate depends smoothly on the other two 
variables (which are angles), then the attractor lies on a 
two-torus, and the considerations of Sec. II apply directly. 
On the other hand, the existence of such a smooth two- 
torus is in question, and this is the main issue addressed 
in Sec. III. A. For the map of Sec. III.B we are able to gen- 
eralize the rigorous approach of Sec. II. A to show that 
for this class of maps Dq = N + 1 while Di — N. In ad- 
dition, numerical experiments are performed to test the 
convergence of dimension computations to these values. 

II. TWO DIMENSIONAL MAP ON A TORUS 
A. Existence of SNA 

We investigate the simplest case of (^ where N — 1 
(5)(«) _> Q'j and the state variable x is one-dimensional. 
Specifically, we take x to be an angle variable tp, so that 
the map operates on a two-dimensional 9-ip torus. Within 
this class we restrict consideration to maps of the form 

0,1+1 = [0,1 + uj] mod 2n, (5a) 
fn+i = [On + Pn + vPiOn, fn)] mod 27r, (5b) 

where uj — 7r(V5 — 1), and P{9,ip) is continuous, dif- 
ferentiable, and 2tt periodic in both of its arguments {9 
and ip). When rj is small enough {\r]\ < 77^), this map 
is invertible. That is, the map is solvable for {9m Pn) 
when {9n+i, Pn+i) is given. We choose a simple func- 
tion P{9, if) — simp for our numerical work. In this case, 
the system is invertible if |»7| < 1. Furthermore, since 
the map is invariant under the change of 77 — > —rj and 
</? ^ (p-|-7r, it is sufficient to consider only the case rj > 0. 

Figure |^ illustrates how a curve C on the 9-p toroidal 
surface is mapped to a curve C" by the map (|^). Note 
that the torus is unrolled in the 9 direction to visualize 
the whole curve C in a two-dimensional plain, but still 
rolled in the tp direction. The curve C circles around the 
torus in the 9 direction, but does not wrap around the 
torus in the (p direction. After one iterate of (||), the curve 
C is mapped to a curve C that wraps once around the 
torus in the (p direction. This behavior comes about due 
to the term 9n on the right-hand side of (^), because 
9 + p + riP{9, p) increases by 2tt as 9 increases by 27r. 
Similarly, applying the map to C produces a curve with 
two wraps around the torus in the p direction, and so on. 

The main results of our numerical experiments and 
rigorous analysis of (^) with \ri\ < rjc are as follows: 

(i) The map (||) has a single attractor. 

(ii) For typical P{9,ip), the attractor has a Lyapunov 
exponent h^p that is negative for rj ^ 0. 

(iii) The attractor has information dimension one for 

(iv) The attractor is the entire 9-p torus and, hence, 
has box-counting dimension two p4| . 



3 



(v) These results are stable to perturbations of the 
system §§0. 

We first establish (ii) using an approximate formula 
for for small 77. Our evidence for (ii) is strong but 
a rigorous mathematical proof is lacking. If we adopt 

(ii) as a hypothesis, then all the other results rigorously 
follow. 

Lyapunov Exponent: A trajectory of the map (|^) 

has two Lyapunov exponents he and hip, where hg = is 
associated with (||a) and /i,^ is associated with (||b). The 
latter exponent is given by the formula, 

h<p^ J H^ + vP<p{d,v)W, (6) 

where = dP/dip, and fi denotes the measure gener- 
ated by the orbit from a given initial point (^oi fo)- 

lih,p > for a particular trajectory, then, since kg = 0, 
the map exponentially expands areas near the trajectory 
in the limit n — > cx). Since the 9-ip torus has finite area, if 
the map is invertible, then there cannot be a set of initial 
points of nonzero area (positive Lebesgue measure) for 
which h;p > 0, and the map thus does not have a chaotic 
attractor. Thus h^p < for typical orbits. 

Furthermore, we argue that h^p < for small nonzero 
rj. We consider first the case 77 = 0, for which (||b) 
becomes (pn+i — (On + ^n) mod 27r. If we initial- 
ize a uniform distribution of orbit points in the 0-ip 
torus, then, on one application of the rj = map, the 
distribution remains uniform. Furthermore, this uni- 
form distribution is generated by the orbit from any 
initial condition. To verify this, we note that the ex- 
plicit form of an 77 = orbit, 0„ — {Oq + nu) mod 27r, 
fn = [</2o + -\- ^{n^ — n)iij] mod 27r, which is shown 
to generated a uniform density in Ref . . We can ob- 
tain an approximation to for nonzero but small 77 by 
expanding ln(l -I- TyP^) in (H) to order 77^ and assuming 
that, to this order, the deviation of the measure fj, from 
uniformity is not significant [d^ sa dO dip /(2Tr)'^]. Using 
ln(l + r]Pip) = rjPip - {l/2)r]^P^ + 0{r]^), this gives 

K ^~lv'<P'> +o{Tf ), (7) 

which is negative for small enough f] 0. Here < P^ > 
denotes the d-ip average of P^, and the order rj term 

is absent by virtue of J^^ Ptpdcp = 0. Since we cannot 
show convergence of an expansion in 77, our result (^ is 
formal rather than rigorous. However, numerical results 
strongly support (0) . Figure || shows a plot of h^p versus 
■q for P{9,ip) — sniLp. Remarkably, Eq. (^ (the straight 
line) describes the numerical data to better than 8 % 
even for 77 as large as 0.5. 

Dimensions of the SNA: For our map the infor- 
mation dimension cannot be less than one due to the 
quasiperiodic 9 dynamics. In addition, the Lyapunov di- 
mension is an upper bound of information dimension [ p^ . 
Therefore, if we accept (ii), h^p < 0, then hg — implies 

(iii) . 



Results (iii) and (iv) quantify the strangeness of the 
attractor. In particular, since the information dimension 
of the attractor is one, orbits spend most of their time on 
a curve-like set; yet, since the box-counting dimension is 
two, if one waits long enough, a typical orbit eventually 
visits any neighborhood on the 9-(p torus. One can get 
a sense of this result from the numerical orbit shown in 
Fig. H, in which a trajectory of length 10^ appears to be 
concentrated along one-dimensional strands [Fig. ||(a)], 
but for the same parameters a trajectory of length 10^ 
fills much more of the 9-ip torus [Fig. ||(b)]. 

We show in Fig. 0(a) a plot of log2 N{e) versus 
log2(l/e), and in Fig. ^(b) a plot of log2(l/_Pi) versus 
log2(l/e). Here N{e) is the number of £ x e boxes (in 9-ip 
space) needed to cover the points from an orbit of length 
T, and pi is the fraction of those orbit points in the ith 
£ X £ box. According to our previous arguments on di- 
mensions, in the limit T — > 00, the points in Fig. ^(a) and 
Fig. ^(b) should follow a straight line of slope two and 
one for small £, corresponding to a box-counting dimen- 
sion of two and an information dimension of one. As is 
commonly found, the box-counting dimension computa- 
tion converge rather slowly with increasing orbit length 
T. Thus, we show plots in Fig. ^ for several different T. 
As can be seen in Fig. ^a), the e range consistent with 
a slope of two (the straight line in the figure) steadily in- 
creases toward smaller £ [larger log(l/£)] as T increases. 
This is contrast with Fig. ^(b), which appears to reach 
a form independent of T that is consistent with a small 
£ slope of one. While the convergence in Fig. ^(a) and 
^(b) is consistent with box-counting and information di- 
mensions of two and one, the slowness of the convergence 
also indicates that a purely numerical determination of 
the dimension values is suspect. 

Topological Transitivity: To establish results (i) 
and (iv), that the attractor of the map is in the whole 
9-p torus, we prove that the map is topologically transi- 
tive: For every pair of open disks A and B, there is a 
trajectory that starts in A and passes through B. This 
property is known to imply that a dense set of initial 
conditions yields trajectories each of which is dense in 
the torus [13. In particular, any attractor, having an 
open basin of attraction, must contain a dense orbit, and, 
hence, must be the entire torus. 

We will show in fact that for every pair of line seg- 
ments Sa — {{9,p) : 9 Ra and (p = pa} and Sb = 
{{9, ip) : 9 e Rb and cp = ipb}, where Ra = {Oa, 9a + Sa) 
and Rb = {9b, 9b + Sb), there is a finite trajectory of M 
that begins on the first segment and ends on the sec- 
ond. (Choosing Sa to lie in A and Sb to lie in B, this 
implies topological transitivity.) In other word, we will 
show that the ?T.th iterate of Sa intersects Sb for some 
positive integer n; see Fig. ^(a). Our strategy is to it- 
erate Sa forward until the union of its iterates covers all 
values of 9 at least once; the number of iterates needed 
is finite and depends only on da- By selecting pieces of 
some of these iterates that cover each value of 9 exactly 
once, we form the graph ip = ga{9) of a piecewise con- 
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tinuous function ga', see Fig. ^(b). Similarly we form a 
graph ip = gb{0) from pieces of backward iterates of 3^. 
Finally, we show that some forward iterate of the graph 
of Qa must intersect the graph of gf,. 

The following is a formal definition of ga- Let Me be 
the map (^). For each 9, let k{9) be the smallest non- 
negative integer for which 9 £ Mg{Ra)- (In Fig. ^(b), 
k{9) = 0, 1, or 2 for all 9.) Let ga{(^) be the (^-coordinate 
of the k{9)-th iterate under M of {M'^^'^'^ ,Lp) E Sa- Then 
the graph ip ~ ga{9) has a finite number of da of disconti- 
nuities. Each contiguous piece of this graph is a forward 
iterate of some piece of Sa- 

Now form the curve Ga by taking the graph of ga and 
adding line segments in the ip direction at each value of 9 
where ga is discontinuous. (We take these segments to lie 
in < (/9 < 27r.) Thus we make Ga a contiguous curve. 
See Fig. |5|(b), which illustrates this construction on for 
a case where da = 3. Notice that, for each n, the nth 
iterate of Ga is also a contiguous curve that consists of the 
graph of a function with da discontinuities, together with 
da "connecting segments" . Define gb and Gb similarly to 
ga and Ga, but in terms of the backward (not forward) 
iterates of the Sb- Let db be the number of discontinuities 
of gb- 

Our goal is to show that for n sufficiently large, the 
nth iterate of Ga intersects Gb for at least da + db + 1 
different values of 9. Then since there are at most da + 
db values of 9 at which one of these two curves has a 
connecting segment, there will be at least one intersection 
point between the nth iterate of the graph of ga and the 
graph of gb- Since the graph of ga consists of forward 
iterates of Sa and the graph of gb consists of backward 
iterates of Sb, some forward iterated Sa will intersect Sb, 
as we claimed. 

Given a contiguous curve C that, like Ga and Gb, is 
the graph of a function of 9 that is continuous except for 
a finite number of values at which C has a connecting 
segment, observe that its image under Alg is a curve of 
the same type (in particular, since the map is one-to- 
one, the heights of the vertical segments remain less than 
27r). Furthermore, because of the 9n term in the ip map 
(||b), the image of C "wraps around" the torus in the ip 
direction one more time than G does as one goes around 
the torus one time in the 9 direction (see Fig. p. 

To formulate what we mean by "wrapping around" , 
define the winding number of G as follows. As 9 increases 
from to 2iT, count the number of times C crosses ip = 
in the upward and downward directions. The difference 
between the number of upward and downward crossings 
is the winding number of C. (The numbers of upward 
and downward crossings may depend on the arbitrary 
choice of as the p value at which to count crossings, 
but their difference does not.) For example, in Fig. ^(b) 
the winding number of Ga is 0. 

Now if two curves Ci and C2, as described above, have 
different winding numbers wi and 'W2, then Ci and C2 
must intersect at least jiui — u'2| times. Because of the 
periodicity of P{9,ip), the winding number of a curve 



must increase by 1 each time the map M is applied. Thus 
for n sufficiently large, the winding number of the nth 
iterate of Ga differs from the winding number of Gb by 
at least da + db + l- Hence the nth iterate of Ga intersects 
Gb for at least da + db + l different values of 9 as desired. 
This establishes claims (i) and (iv). 



Notice that the argument above does not depend on 
the specific form of P{9,(j)), only that it is continuous 
and periodic and that rj is sufficiently small {\r]\ < rjc) 
that the map (||) is one-to-one. This independence of the 
results from the specific form of P{9, cj)) implies that the 
results are stable to system changes [our claim (v)] that 
preserve a quasiperiodic driving component (Sa). 



Discussion: The possible existence of SNA's was orig- 
inally pointed out in and many numerical explo- 
rations of the dynamics on attractors that are appar- 
ently strange and nonchaotic have appeared. Recently, 
there has also been rigorous results on the mathematical 
properties that SNA's must have if they exist In 
spite of these works, a very basic question has remained 
unanswered: Can it be rigorously established that SNA 's 
generically exist in typical quasiperiodically forced sys- 
tems? This is an important issue, because, although the 
numerical evidence for SNA's is very strong, perhaps the 
attractors observed are nonstrange with very fine scale 
structure (rather than the infinitesimally fine scale struc- 
ture of a truly strange attractor). Also, there might 
be the worry that the numerical evidence is somehow 
an artifact of computational error. Our proof of topo- 
logical transitivity, combined with the hypothesis that 
< 0, answers the question of the typical existence 
of SNA's (affirmatively) for the first time ([|l3| contains 
a preliminary report of our work). The only previous 
work rigorously establishing the existence of a SNA is 
that appearing in the original publication on SNA's 0] 
and in |21[| . These proofs, however, are for a very spe- 
cial class of quasiperiodically forced system such that 
an arbitrarily small typical change of the system puts 
it out of the class. Thus this proof does not establish 
that SNA's exist in typical quasiperiodically forced sit- 
uations. In order to see that nature of this situation 
with respect to Refs. |^ and we recall the example 
treated in Ref. In that reference the map consid- 

ered was Xn+i = 2A(tanha;„) cos6'„ = /(x„, 0„), with 9n 
evolved as in Eq. (js^) . It was proven in Q that this map 
has a SNA for A > 1. However, the map has an invariant 
set, namely, the line x = 0, 9 in [0, 27r), and this fact is 
essential in the proof of Ref. On the other hand, 

the existence of this invariant set does not persist under 
perturbations of the map. Thus, if we perturb f{x, 9) 
to f(x,9) + eg(x,9), the invariant set is destroyed, even 
for small e, for any typical function g{x, 9) (in particu- 
lar, an arbitrarily chosen g{x,9) is not expected to satisfy 
5(0,0) =0). 
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B. Origin of SNA's: Rational Approximation 

Using rational approximations (RA's) to the quasiperi- 
odic forcing, we now investigate the origin for the ap- 
pearance of SNA's in (||) for P{9,(p) = sin ip and uj = 
7r(\/5 — 1). For the case of the inverse golden mean 
il = w/27r, its rational approximants are given by the ra- 
tios of the Fibonacci numbers, fife — Fk-i/Fk, where the 
sequence of {Fk} satisfies Fk+i = Fk + Fk-i with Fq ~ 
and Fi = 1. Instead of the quasiperiodically forced sys- 
tem, we study an infinite sequence of periodically forced 
systems with rational driving frequencies ujk- We sup- 
pose that the properties of the original system may be 
obtained by taking the quasiperiodic limit fc — > oo. 

For each RA of level k, a periodically forced map 
with the rational driving frequency Qk has a periodic or 
quasiperiodic attractor that depends on the initial phase 
Oq of the external force. Then we take the union of all 
attractors for different to be the fcth RA to the attrac- 
tor in the quasiperiodically forced system. Furthermore, 
due to the periodicity, it is sufficient to obtain the RA 
by changing 6*0 only in an basic interval S [0, l/i^fc), 
because the RA to the attractor in the remaining range, 
[l/Fk, 1), may be obtained through {Ft — l)-times iter- 
ations of the result in [0,1/ Fk). For a given k we call 
the periodic attractors of period Fk the "main periodic 
component" . The first column of Fig. ^ shows that the 
Lebesgue measure of the main periodic component (de- 
noted by the solid line) becomes dominant as the level k 
increases (i.e., the fraction of the 9 axis corresponding to 
the non-periodic, gray area decreases). By iterating the 
RA in the basic interval of 6, we obtain the RA in the 
whole range of 6, as shown in the second column of Fig. |^. 
As k increases, the whole RA becomes more similar to 
its quasiperiodic limit given in Fig. |^. 

We first note that for 77 = the RA to the regular 
quasiperiodic attractor consists of only the quasiperi- 
odic component. However, as rj becomes positive peri- 
odic components appear via phase-dependent {i.e., Oq- 
dependent) saddle-node bifurcations. As an example, 
we show the 6th RA for rj = 0.3 in Fig. |^. Here the 
quasiperiodic component is plotted in the gray. "Gaps" 
in the gray quasiperiodic regions are occupied by periodic 
attractors. As examples, we explicitly show the main 
period- (^6 — 8) and minor period-SFg components in 
Figs. 0(a) and 0(b), respectively. At both ends of each 
gap, a pair of stable (denoted by a solid line) and unsta- 
ble (denoted by a dashed line) periodic orbits appear via 
a phase-dependent saddle-node bifurcation. Figures ||(a) 
and ^(b) show the saddle-node bifurcation curves in the 
OFk-t] plane, at which the main periodic components with 
period Fk are born. It can be easily seen that for a given 
r] the width of the main gap (occupied by a period- i*!- at- 
tractor) becomes larger as k increases. Quantitatively, it 
is found that the Lebesgue measure /i^ in for the main 
periodic component becomes dominant as k increases; 
i.e., the Lebesgue measure (1 — fik) of the complemen- 
tary set decreases exponentially with Fk; 1 — /ifc ~ e""^*", 



where a = 0.013, as shown in Fig. ||(c) for rj = 0.3. 

In what follows we use the RA's to explain the origin of 
the negative Lyapunov exponent and the strangeness 

of the SNA. For a given level k of the RA, let h^^'^ (9) 
denote the Lyapunov exponent of the attractor corre- 
sponding to a given 0. Thus hlp \e) = for 6* in the 
quasiperiodic range (gray regions of Figs. |^ and 0) and 

/i^*^'' (9) < for 6 in the periodic range (gaps in the gray 
regions). Since the attractor with irrational uj generates 
a uniform density in [see Eq. (|^)] , we take the order-fc 
RA to the Lyapunov exponent h^p to be 

<hf>^^rh'^^\0)d0. (8) 

For 77 > 0, due to the existence of periodic components, 
< h\^^ > is negative. As 77 increases for a given level k, 
the Lebesgue measure in for the periodic components 
increases, and hence hlp"^ (0) becomes negative in a wider 
range in 0, as shown in Figs, ^(a) and||(b) for level fc = 6. 
Thus, as ?7 increases < h''^^ > decreases [see Fig. ^(c)]. 

In addition, we note that as the level k increases, the 
RA to the Lyapunov exponent h^p converges rapidly to 
its quasiperiodic limit [represented by the solid line in 
Fig. H(c)]. For comparison, the approximate analytic re- 
sult for (i.e., h^p = —77^/4) is also given (see the dashed 
line). Consequently, for any nonzero rj the attractor in 
the quasiperiodic limit has a Lyapunov exponent whose 
value decreases as rj is increased. 

We now discuss the strangeness of the attractor in the 
quasiperiodic limit for 77 = 0.3. In the quasiperiodic case, 
we have seen (Fig. ||) that a typical trajectory seems to 
fill the whole torus densely, but, unlike the case of the 
regular quasiperiodic attractor, it appears to spend most 
of its time on a set of ID strands. We identify these 
apparent ID strands with the fc — > 00 limit of the main 
periodic component. Although, as k becomes larger, the 
Lebesgue measure of the quasiperiodic region approaches 
zero [Fig. ^(c)], these quasiperiodic region become dense 
in 0. Since each quasiperiodic region fully covers the ip 
interval [0, 27r), the attractor is expected to occupy the 
entire 9-ip torus, and, hence, it is expected to have a 
box-counting dimension of 2. 



III. HIGH DIMENSIONAL MAPS 

A. Radial Perturbations of the Torus Map 

We now show that stability to perturbations applies in 
addition if the system is higher dimensional. In partic- 
ular, we discuss the case of a three-dimensional system 
with an attracting invariant torus, and allow perturba- 
tions of the toroidal surface. Consider the following map 
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on R^: 

Sn+i = [fi'n + ^] mod 27r, (9a) 

•fn+i = [Qn + ^n+riP{Qn,^n,rn)] mod 27r, (9b) 

r„+i = Ar„ + pQ{9n, fn)- (9c) 

Here and are coordinates on a torus embedded in R'^, 
as m Fig. [l[ and r is a coordinate transverse to the torus, 
with r = representing the unperturbed {\ = p ~ 0) 
torus. The parameters uj and 77, and the dependence 
of P on and (p, have the same properties as for map 
(^, and Q is continuously differentiable and 27r periodic 
in 9 and (p. When A and p are small, Eq. (^ maps a 
neighborhood of the torus r = into itself, and when 
p = the torus r = is invariant and attracting. It then 
follows from classical results on the perturbation of in- 
variant manifolds |^ that, for A and p sufficiently small, 
the map (^ has a smooth attracting invariant manifold 
r = /(0, (/?) near the torus r = 0. On this attractor, 
the map (g) reduces to a map of the form (||), with 
P{0,ip) = P[9,(p,f{e,ip)]. Thus statements (i)-(v) above 
apply also to the attractor of the three-dimensional map 




The above arguments depend on the existence of a 
smooth invariant torus on which the attractor is located, 
and this is guaranteed if A and p are sufficiently small. 
We now show numerical evidence for the existence of a 
smooth invariant torus for values of A and p that are ap- 
preciable. We consider the example P{0, ip, ?') = sin (p 
and Q{9,(p,r) = sin(r + cp). We numerically obtain 
two-dimensional plots of intersections of the invariant 
torus with the surfaces ip — n [Fig. |l^(a)] and If ^ TT 
[Fig. |l^(b)]. First consider the case tp = TT. Our nu- 
merical technique is as follows. We choose an initial 
value (6*0,^0 = tt) and obtain (6'„„, (/3_„) by iterating 
(^) and (|yb) backward n steps. Since /i^ ~ In A < 
(when p <C A < 1), r_„ — s- ±00 if tq is not on the torus. 
In other words, if r_„ = 0, then tq is on the torus. Thus, 
we choose r_„ = and iterate (r_„, 0_„, c/j_„) forward 
n steps to {rQ,dQ,ipQ = tt). By varying 0o, we obtain 
the graph, to^Oq), of the torus intersection with tp = n. 
Similarly, choosing (^o — ""i Vo) a-nd iterating the map 
backward, and then forward, we can obtain ro^tpo) 
of the torus intersection with 9 = tt. (For our numeri- 
cal experiments, we set n — 25.) As shown in Fig. 
a smooth invariant torus exists in the parameter region 
where p and A are appreciable. 

In Figs. pl](a) and (b) we show dimension computations 
for the map (^) with P{9, tp, r) — sin((p) and Q{9, ip, r) = 
sin(r -|- <p). [In this three dimensional case we employ e 
edge length cubes in (0, tp, r)-space.] As for Fig. |[ we 
observe from Fig. |l^, that the results are consistent with 
slow convergence to the predicted dimension values of 2 
and 1 as the orbit length T is increased. 



B. Map on a High Dimensional Torus 

In Section II. A we proved that (||) is topologically tran- 
sitive. Here we show how this argument can be mod- 
ified to higher dimensional maps that include TV > 1 
quasiperiodic driving variables 6'^-'^', 9'^'^\ • • • , 6'^^^ For 
exposition we assume N — 2, but the argument is virtu- 
ally identical for all N. 

In particular, we consider a map of the form, 

= [0(^1) + mod 2^, (10a) 

9'^l^ = [0(^2) ^ ^(2)j j^qJ 2^, (10b) 
Vn+l^[q^'^9iP+q('^9l^^+^,,+ (lOc) 
r,P{9^^\9l^\ipn)] mod 271, 

where cj^^^ and w*^^) are incommensurate, {q^^\q^'^'^) is a 
pair of integers different from (0,0), and P {9'^'^\ 9'''^\ >p) 
is continuous, differentiable, and 2tt periodic in all of its 
arguments (61(1), 61(2), and p). We assume without loss of 
generality that q^^) ^ 0. 

Let Ra = {(0(1), 0(2)) . oil) ^ ^(1) ^ (g(i) ^ 

6a) and 0(2) < 0(2) < (0(2) +,5,)} and Rt - {(0(i),0(2)) : 
0^'^ < 0(1) < (0^') + Sb) and 0f ) < 0(2) < (0f ) + 5b)} 
be two arbitrary squares in the 0(i)-0(2) torus, and let 
Sa = {(0(1), 0(2), V?) : (0(1), 0(2)) e Ra and ip = p^} and 
Sb = {(0(1), 0(2), (p) : (0(1), 0(2)) g Rb and ip = ipb} be a 
pair of square segments, where tpa and tpb are arbitrary. 
As before, we will show that there is a finite trajectory 
that begins on Sa and ends on Sb- 

In this case, we proceed by iterating Ra forward until 
the union of its iterates cover all points (0(i) , tt) at least 
once (see Fig. |l2|). The number of iterates needed is 
finite. Then we select pieces of these iterates that single- 
cover a thin strip Da = {(0(i), 0(2)) : tt < 0(2) < tt + Sa} 
with rectangles of width Sa- From the corresponding 
pieces of the corresponding iterates of Sa, we form the 
graph (p = ga{9'^^\9'^^^) of a piecewise continuous func- 
tion ga defined on Da- Similarly we form a graph 
ip = gb{9^^\9^^^) on a strip Db from pieces of backward 
iterates of Rb x {ipb}- As before, we will show that some 
forward iterate of the graph of must intersect the graph 
of gb. 

Next, form the strip Ga by taking the graph of ga and 
adding "connecting faces" at each of the da values of 0(i) 
where ga is discontinuous, so as to make Ga a contiguous 
strip. The construction of Ga is essentially as shown in 
Fig. ^(b), except that it now has some thickness in the 
0(2) direction (not shown). For each n, the nth iterate of 
Ga is also a contiguous strip that consists of the graph 
of a function with da discontinuities in the 0(i) direction, 
together with da connecting faces, over a strip in the 
0(i)-0(2) torus of width Sa in the 0(2) direction. Notice 
though that the strip moves a distance Li;(2) in the 0(2) 
direction with each iteration. Define gb and Gb similarly 
to ga and Ga, but in terms of the backward iterates of 
Sb, and let db be the number of values of 0(i) at which 
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gb is discontinuous. 

As before, we can define the winding number of strips 
like Ga and Gb, representing the net number of times the 
strip wraps in the f direction as 9^^^ increases from 
to 2tt. The winding number can be computed for any 
fixed value of 0^^^ and does not depend on that value. 
With each iteration of (|lO|), the winding number of such 
a strip changes by q^^^ ^ 0. Therefore for n sufficiently 
large, the winding number of the nth iterate of Ga differs 
from the winding number of Gb by at least da + db + 1- 
Furthermore, by increasing n if necessary, we can ensure 
that the domains of these two strips intersect; that is, 
they have a common value of 0'^^^ . Then for that value of 
O^^^f , it follows as before that the nth iterate of the graph 
of ga (without the da connecting faces of Ga) and the 
graph of gb (without the db connecting faces of Gb) must 
intersect as claimed. 

Our numerical experiments also give a sense of the 
above proof. In order to obtain two-dimensional plots, 
we count points when the trajectory passes through a 
thin slab of width 5 <C 1 containing the 6'^*'-iy9 surface. 
(For our experiments, the width of the slab is 5 = 0.01.) 
Figures ^ show numerical approximations, thus, ob- 
tained, to the intersections of the attractors with the 
6»(2) = TT surface [Fig. [l3|(a)] and the 6l(i' = n surface 
[Fig. |l3|(b)]. Both figures look similar to the figure for 
two dimensional case (Fig. |^) . We also show in Fi g. H (a) 
a plot of log2A^(e) versus log2(l/e), and in Fig. |14|(b) a 
plot of 'Spi log2(l/pi) versus log2(l/£). According to the 
above proof, the box-counting dimension should be three 
and the information dimension is two. However, due to 



the slowness of the convergence of slopes, the dimension 
values are not clearly evident, although, as in the two- 
dimensional case (Sec. II. A) the results are suggestive. 



IV. CONCLUSION 

In this paper we addressed the existence of robust 
strange nonchaotic attractors. In particular, we pro- 
vided rigorous analysis for the two-dimensional map (|^) 
in Sec. II. A and for the {N + l)-dimcnsional maps of the 
form of (H) in Sec. III.B. In addition, we have used a ratio- 
nal approximation technique (Sec. II. B) to investigate the 
dynamical origin of SNA's, as well as to gain additional 
understanding on why Di = 1 and Dq = 2. In Sec. III. A, 
we show that the stability to perturbations (robustness) 
continues to apply in systems [e.g., (^] where there can 
be an attracting torus. We also carried out calculations 
to see how our rigorous dimension results are manifested 
numerically. Our results confirm the existence of SNA's 
as a generic phenomenon of quasiperiodically forced sys- 
tems. 
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FIG. 1: Torus unwrapped in the 9 direction = and 9 = 2n 
are identified with each other). The map (p|) takes the curve 
C to the curve C'. 



FIG. 2: Lyapunov exponent versus r;^ . For each t], the data 
plotted as open circles were computed from 10^ iterations of 
the map with ui = tt {V5-1) and P{9,ip) — sm ip. 

FIG. 3; Trajectory of the map (P) with lu = n{\/5 — 1), tj = 
0.3, and P{9, ip) = sin (p. In each case 9o = tpo — and lO** 
points of the trajectory are computed before plotting; in (a) 
the next 10'* points are plotted, while (b) shows 10^ points. 



FIG. 4: Dimension computations for (|5|) with 7] = 0.3, 

UJ = TV (75 - 1), and P{9,ifi) — sin (p. In (a) the dashed 
line has slope two, while in (b) it has slope one. In each 
graph, the curves from lowest to highest represent T = 
10^, lO**, ■ ■ ■ ,10^"; in (b) the final five curves are virtually 
identical. 



FIG. 5: (a) The nth iterate of Sa intersects Si,. The nth 
pre-iterate of this intersection point (denoted p) is a point 
on Sa that goes to Sb in n iterates, (b) Sa plus its first two 
iterates, M{Sa) and M^{Sa), cover the entire 2% range of 9. 
M{Sa) and M^{Sa) are shown as thin lines. The curve Ga, 
which includes Sa, pieces of A{{Sa) and M'^{Sa), and vertical 
segments connecting these pieces, is shown as a dark thick 
line. 

FIG. 6: RA's for 77 = 0.3. The levels are it = 6 in (a) and (b), 
fc = 8 in (c) and (d), and fc = 11 in (e) and (f). In the first 
column the RA in the basic interval of 9 is given, while in the 
second column the RA in the whole range of 9 is given. The 
quasiperiodic component is represented in gray dots and the 
main periodic component is denoted by the solid line. 

FIG. 7: 6th RA for — 0.3. The quasiperiodic component is 
shown in gray, and the stable and unstable periodic orbits in 
the gaps are denoted by solid and dashed lines, respectively. 
Main period- and minor period-3-F6 components are shown 
explicitly in (a) and (b), respectively. 



FIG. 8: Phase-dependent saddle-node bifurcation lines for the 
main periodic components. The cases of the level fc = 6, 9, 12 
are shown in (a), and other cases with fc = 7, 8, 10 are given 
in (b). (c) Plot of ln(l — jj,k) vs. Fk for rj = 0.3. Solid points 
denote the data for level fc = 6, . . . , 12, which are well fitted 
with a dashed straight line with slope a = 0.013. 

FIG. 9: Plot of h':f\e) vs. OFe for (a) -q = 0.1 and (b) r] = 0.3. 
(c) Plot of < /i^''' > vs. rj for the three levels fc = 6, 7, and 8. 
The solid and dashed lines denote the Lyapunov exponents 
in the quasiperiodic limit that are obtained numerically and 
analytically, respectively. 
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FIG. 10: Attractors in two-dimensional surfaces, (a) r versus 
6 &t if = n surface and (b) r versus ip a,t 6 — n with p — 0.5, 
A = 0.5, and r) = 0.3 {h^ = -0.024 and hr = -1.370). 



FIG. 11: Dimension computations for with rj = 0.3, A — 
0.5, p = 0.5, io — ■k{\/5 — 1), P{6, ip) = sin (p, and Q{0, tp, r) — 
sin(r + v')- In the dashed hne has slope two, while in (b) 
it has slope one. In each graph, the curves from lowest to 
highest represent T = 10^, 10*, • ■ • , 10®; in (b) the final four 
curves are virtually identical. 



FIG. 12: Construction of the domain Da (shaded region) of 

9a- 



FIG. 13: Trajectory of the map (|l|) with oj*^' = 7r(^- 1), 
J^'^ = 2^l{^/2 - 1), r? = 0.3, g'^' = = 1, and 

P(6l(^>,6l(2),(^) = sin(6l(^) +ip). In each case 0*^' = 6^^ ^ = 
<Po = and 10* points of the trajectory are computed before 
plotting the next 10^ points in the slice of (a) tt- 5/2 < 0^^^ < 
n + (5/2 and (b) tt - 5/2 < 9^^'> < n + 5/2, where S = 0.01. 



FIG. 14: Dimension computations for ( |lC| ) with r; = 0.3, 
cjti) = 7r(75 - 1), = 27r(V2 - 1), and P{e^^\6^^\p) = 
sin(6''-^-' + ip). In (a) the dashed line has slope three, while in 
(b) it has slope two. In each graph, the curves from lowest to 
highest represent T = 10*, 10^ ■ ■ • ,10®. 
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